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Synchrotron radiation of dissipative solitons in optical fiber cavities 
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New resonant emission of dispersive waves by oscillating solitary structures is considered analyt¬ 
ically and numerically.The resonance condition for the radiation is derived and it is demonstrated 
that the predicted resonances match the spectral lines observed in numerical simulations perfectly. 
The complex recoil of the radiation on the soliton dynamics is discussed. 
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Since the pioneering work by Akhmediev and Karlsson 
[I| the resonant radiation of optical solitons propagating 
in fibers with high order dispersion has been attracting 
much of attention because of both the fundamental in¬ 
terest and practical importance HQ. The condition of 
the resonance is the equality of the soliton velocity to the 
phase velocity of a dispersive wave and in this way the 
radiation is similar to the effect of Cherenkov radiation 
of charged particles. A similar effect of transitional ra¬ 
diation of Bragg solitons was later discovered in periodic 
media Q- 

The aim of the present Letter is to address the reso¬ 
nant radiation of oscillating optical solitons and to reveal 
the analogy of this radiation to the synchrotron radia¬ 
tion of charges appearing when the velocity of relativistic 
charges varies periodically in time [5]. 

It should be noted here that from the point of view 
of the field theory the synchrotron radiation can be in¬ 
terpreted as an excitation of resonant waves by an oscil¬ 
lating localized source moving at a relativistic velocity. 
The periodical variations of the velocity of a relativistic 
charge is just a way to obtain the right driving force in 
the equation for the electromagnetic waves. 

Aiming to develop the analogy between synchrotron ra¬ 
diations of the charges and solitons it is much simpler and 
more instructive to consider one-dimensional oscillating 
solitons rather than two-dimensional nonlinear localized 
structures relativistically moving along one direction and 
wobbling in the transverse direction (though the latter is 
also possible). 

The effects of Cherenkov and synchrotron radiation 
are widely used in vacuum electronics, in particular in 
forward- and backward waves oscillators. The optical 
counterparts of Cherenkov effects can play an important 
role in the generation of new optical frequencies by the 
solitons propagating in fibers with high order dispersion. 
Indeed, since 1995 when Cherenkov radiation of optical 
solitons was reported for the first time this effect has 
been actively studied both experimentally and theoreti¬ 
cally, especially in the context of the generation of optical 
supercontinuum HQ. Recently a number of papers were 
published on Cherenkov emission shed by shock waves jb] 
and ultra-short pulses HQ- 

All these effects are important for the explanation 
of the dynamics of solitons and domain walls in wave- 
guiding systems. The effect of synchrotron emission of 


oscillating solitons can also affect the dynamics of the 
solitons strongly and thus is of interest from both funda¬ 
mental and practical points of view. 

Optical solitons in fibers seem to be the best candidates 
for experimental observation of the discussed effect of 
synchrotron emission, however similar effects can take 
place in polariton condensates, plasmas, hydrodynamics 
and other systems of different physical origins. 

The resonant emission is possible not only in the con¬ 
servative but also in weakly dissipative systems. In this 
Letter we consider the radiation in annular dissipative 
optical cavities with external pump. The advantage of 
dissipative systems is that the resonant radiation can oc¬ 
cur there in a stationary regime when it is not disguised 
by any transitional processes. 

Our choice is also motivated by the fact that these 
systems are not only a convenient test bench for inves¬ 
tigation of synchrotron radiation but are of interest on 
their own and have been actively inv estigate d in recent 
years theoretically and experimentally [12l - [24{ . 

Propagation of oscillating solitons in the presence of 
high order dispersion was studied experimentally and nu¬ 
merically and very reach dynamics of the solitons was 
reported [25]. Cherenkov radiation of oscillating cavity 
solitons was recently considered in [26|] but the effect of 
the synchrotron radiation was out of the scope of that 
paper. In the present Letter we focus on synchrotron 
radiations specific for oscillating solitons. 

To describe the cavity we adopt a well known Lugiato- 
Lefever model in the form as it is written in [23| 

d t A = -(1 + iO)A + i\A\ 2 A + id 2 x A + d s d*A + P (1) 

where 0 is the detuning between the cavity resonance 
and the frequency of the pump, P is the amplitude of 
the pump and d% is the third order dispersion. To make 
the resonance spectral lines well pronounced we choose 
the detuning to be equal to a realistic value 6 = 15 which 
is just slightly greater then the values used in (23|. The 
boundary conditions for annular cavities are periodic but 
it is checked that for the for large cavities considered in 
the present Letter the boundary conditions do not affect 
the dynamics of the solitons. 

Equation m has bright soliton solutions that loose sta¬ 
bility through Andronov-Hopf bifurcation leading to the 
formation of an oscillating soliton mm for the pumps 
exceeding a threshold value. For the chosen parame- 
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FIG. 1: (Color online) Panel (a) shows the temporal evolution 
of the field of an oscillating soliton for d 3 = 0.02 and P — 8 . 
Typical spectra of the fields are shown in panel (b) for d 3 = 
0.02, d 3 = 0.04 and d 3 = 0.06 by the black, blue and red lines. 
The synchrotron radiation spectrum for d 3 = 0.02 is shown 
in the inset. 


ters the threshold pump is P t h ~ 7. Aiming to consider 
the behaviour of oscillating solitons we choose the pump 
above the threshold P = 8 and consider dynamics of the 
soliton. 

The stationary temporal evolution of the intensity of 
the field is shown in panel (a) of Fig.[l]for relatively small 
value of third order dispersion d 3 = 0 . 02 . It is clearly seen 
that the soliton oscillates in time. The typical spectra of 
the field are shown in panel (b). The resonant emission is 
very weak for d 3 = 0.02 but still clearh^seen at k « 50 as 
predicted by the resonance condition [26j . Much stronger 
radiation is seen at k « 28 and k ~ 20 for d 3 = 0.04 and 
d 3 = 0.06. 

A very important fact is that the spectrum of the reso¬ 
nant radiation consists not of a single line as it should be 
in the case of Cherenkov radiation but of a series of well 
resolved peaks. Below we will show that, indeed, one of 
these lines corresponds to Cherenkov radiation whereas 
the other lines can be refereed as synchrotron radiation. 

It is worth mentioning here that narrow radiation lines 
at low k overlapping with the soliton spectrum are also 
seen in panel (b). These lines can be understood as 
non-relativistic radiation of the oscillating solitons and 
thus are analogous to, for example, cyclotron radiation 
of charged particles. 

It is worth mentioning here that the equation (pQ) is of 
course not invariant in respect of Lorentz transformation 
and in this Letter we use the term ” relativistic radiation” 
only in the sense that the velocity of the soliton is close 
to the phase velocity of the dispersive waves. 

To understand the observed phenomenon let us derive 
the conditions of resonant radiation of oscillating soli¬ 
tons. Here we only sketch the derivation, the complete 
analysis will be presented elsewhere. 

Assuming that the radiation is weak it is possible to 
look for a solution in a form A = Aq + u[t,x\ where 
Ao describes the soliton oscillating with period T and 
moving with the velocity v; the term u accounts for the 
small correction. Here and below the arguments of the 
functions are given in square brackets. 


It is convenient to write an equation for u in vector 
form u - (Re u, Im u) T 

d t u + Lu = f. (2) 

The linear operator L has the form 

f _ f 1 + Im Al + d 3 dl 2 1-ReAl-0 + d 2 x \ 
\-21-ReAl + 0-dl 1 - Im A 2 0 + d 3 dl ) ’ 

where I = |An| 2 . The right hand side is given by f = 
(Re /o, Im /o) T with 

fo = P ~ (1 + i0 - i\A 0 \ 2 - idl - d 3 dl)A 0 - d t A 0 . (3) 

The left hand side of equation ([2]) describes propaga¬ 
tion of waves on the background hosting the soliton. In 
the conservative case the resonance of the right hand side 
/ with a delocalized eigenmode of the medium results in 
the formation of a continuous wave propagating away 
from the soliton. In the dissipative case the length of the 
radiation tail generated by the resonance is, of course, 
always finite but can be very large for low losses. 

The oscillations of the soliton make the coefficients in 
L to be periodic functions of time. It means that the 
eigenfunctions of the operator L are Bloch functions in 
time. The important fact, however, is that far from the 
soliton the problem becomes uniform in both space and 
time and thus the asymptotic of the delocalized Bloch 
modes ^ at x —»• ±00 are just plane waves ^ e ^t-tkx 
parametrized by the frequency uj and the wave vector k. 
So the dispersion characteristics if the delocalized eigen- 
modes is nothing else but the dispersion characteristics 
of plane waves uj = c o± [k] on the background 

U± = i + d 3 k 3 r ± vy? + fc 2 -/ o )(0 + fc 2 —3Io), (4) 

where Jo = /[Too] is the intensity of the background. 

To proceed let us represent the right hand side of (|2]) 
in the form / = Ui[x — vt] • e lluJot where uoo = 2tt/T is 
the fundamental frequency of the soliton. Then each of 
the term Ui in the sum can be represented as series over 
the eigenfunctions of the operator L. From the fact that 
all functions Ui[x — vt\ are moving without changing their 
shapes it follows that the frequencies of the driving force 
components depend on their wave vectors as uj = kv-\-luj$. 

The emission of radiation takes place when a harmonic 
of the right hand side / is in resonance with an eigenwave 
of the medium. This gives the condition of the phase 
synchronism 


W±(kr) = k r V + luj 0 - ( 5 ) 

For l = 0 this is the condition of Cherenkov radiation. If 
the condition is satisfied for l 7 ^ 0 then the radiated wave 
is excited by Z-th temporal harmonic of the oscillating 
soliton. 

Let us now discuss how the derived resonance condi¬ 
tions explain the results of the numerical experiments. To 
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FIG. 2: (Color online) Panel (a) shows the spectrum at low 
k of the oscillating (thicker black line) and non-oscillating 
solitons (thinner blue line). The spectrum in the vicinity of 
Cherenkov resonance is shown in panel (b). The graphical 
solutions of the resonance conditions ([5]) are shown in pan¬ 
els (c) and d), the green line corresponds to the Cherenkov 
resonance and the red ones to the resonances with l ^ 0. Res¬ 
onances marked as 1+, 1_ on panel (a) and by 1 on panel (b) 
are given by different crossings of the same red line with the 
dispersion characteristic of the linear excitations on the soli- 
ton background; the same is for the resonances marked by 2. 
The third order dispersion is cfe = 0.02, the pump is P = 8. 

do this we study the spectrum and the resonance condi¬ 
tions at low k and at k close to Cherenkov resonance. We 
start with the resonances at small wave vectors k. The 
spectrum obtained in numerical simulations is shown in 
panel (a) of Fig. [2] The spectrum looks like a bell-shaped 
function with the narrow line at k = 0 corresponding to 
the background. Comparing spectra of non-oscillating 
and oscillating solitons one can notice that in the latter 
case several additional maxima appears on the spectral 
curve. 

The graphical solutions of the resonance condition (j5j) 
are shown in panel (c). The soliton velocity and the 
frequency of the soliton was extracted from numerical 
simulations. It is seen for \l\ >2 the resonance condition 
is met. Comparing the positions of the crossings with 
the positions of the small local maxima of the spectrum 
we see that they match perfectly. The phase velocity of 
these resonant modes are much higher comparing to the 
soliton velocity and thus developing the analogy between 
the radiation of solitons and moving charges we can refer 
these resonances as cyclotron radiation of solitons. 

It is worth noticing that the emitted waves with k < 0 
propagate with the group velocity greater than the veloc¬ 
ity of the soliton and so the radiation appears in front of 
the soliton. The emitted radiation with k > 0 propagates 
in the opposite direction and so the radiation appears be¬ 
hind the soliton. Topologically it is obvious that if there 


is a resonance at a positive k then there is a resonance at 
a negative k and so cyclotron radiation always appears 
simultaneously in front and behind the soliton, see pan¬ 
els (b) and (c) of Fig. [3] showing the field distributions of 
the oscillating solitons. 

Now let us consider the radiation with the frequen¬ 
cies close to the Cherenkov resonance. The crossing of 
the green line (l = 0) with the dispersion characteristics 
in panel (b) gives the position of Cherenkov resonance 
matching one of the spectral lines observed in the numeri¬ 
cal experiment perfectly. The other resonances (/ ^ 0) fit 
the other spectral lines. So we can explain the spectral 
structure of soliton radiation as a number of lines pro¬ 
duced by different temporal harmonics of the oscillating 
soliton. 

It is obvious that the discussed waves have phase ve¬ 
locities close to the velocity of the soliton (for Cherenkov 
resonance they are exactly equal) and thus these radia¬ 
tions can be interpreted as relativistic radiation of the 
oscillating soliton, ie synchrotron radiation. In the con¬ 
sidered case Cherenkov and synchrotron radiations al¬ 
ways have group velocities higher than the velocity of 
the soliton and so the radiation must appear in front of 
the soliton. 

It is instructive to compare the radiation of oscillating 
and not oscillating solitons. The blue curves in panels 
(a) and (b) show the spectrum of a non oscillating soliton 
forming in the case of a weaker pump P = 6. The ve¬ 
locities of the solitons are practically the same for P = 6 
and P = 8 and so the Cherenkov spectral lines of the os¬ 
cillating and non-oscillating solitons practically coincide. 
However it is clearly seen that the non-oscillating soliton 
emits neither cyclotron no synchrotron radiation. 

In the case of third order dispersion the synchrotron 
and Cherenkov resonances always appear together, but 
it is important to remark that there are dispersion char¬ 
acteristics allowing for synchrotron resonances in the ab¬ 
sence of the Cherenkov one. 

To shed more light on the resonant emission we con¬ 
sider the shapes of the radiation fields. Panel (a) of Fig. [3] 
shows stationary Cherenkov radiation of non-oscillating 
soliton pumped by P = 6 when third order dispersion 
is equal to d% = 0.02. The radiation contains only one 
mode, see the spectrum in panel (b) of Fig.[2j and decays 
exponentially with the rate l/\v g — v\. 

For stronger pump P = 8 the soliton develops oscil¬ 
lations and, as discussed above, the cyclotron radiation 
appears in front and behind the soliton. This radiation 
is clearly seen in panel (b). The radiations having small 
group velocity decays along x quickly, so that this radia¬ 
tion is seen only in the vicinity of the soliton. 

Cherenkov and synchrotron radiations with large k 
have bigger group velocities and so decay slower then 
cyclotron radiation. This radiation marked by SR is 
weak and barely visible for ds = 0.02 but becomes much 
stronger for d% = 0.06, see panels (b) and (c) 

A very important fact is that the radiation of an os¬ 
cillating soliton is pulsating whereas the radiation of a 
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FIG. 3: (Color online) Radiation field of the non-oscillating 
soliton is shown in panel (a), P = 6 cfa = 0.02. The radi¬ 
ations of the oscillating solitons (P gg 8 ) are illustrated in 
panels (b) and (c) for d 3 = 0.02 and d 3 = 0.06. The syn¬ 
chrotron and cyclotron radiations are marked by ”SR” and 
” CR” correspondingly. 



FIG. 4: (Color online) Black and green lines in panel (a) 
show the dynamics of the solitons in dtlmax-Imax plane for 
cfe = 0.06 and cfa = 0.02 correspondingly. The evolution of 
the spectrum of an oscillating soliton is shown in panel (b) 
for g ?3 = 0.06, panel (c) shows the spectra of synchrotron 
radiation at two different times. The pump is P — 8 . 


non-oscillating soliton is a monotonically decaying con¬ 
tinuous wave. The pulsations can be explained by the in¬ 
terference of synchrotron emissions at close frequencies. 
The distances between the synchrotron resonances can 
be estimated as A k = uj o/\v gs — v\ where v gs = dkuj(k r ) 
is the group velocity of the radiation. Then the length of 
the radiation pulses is 2nv gs / Auj giving a good estimate 
of w 15 for ds - 0.06. 


Now we turn our attention to a recoil from the ra¬ 
diation on the oscillating solitons. From numerics we 
extract the dependence of maximum intensity of the soli¬ 
ton Imax [t\ on time and calculate its temporal derivative 
dtlmax[t \• These dependencies parametrically define the 
curves in I max [t]-d,L max [£] plane, see (a) of Fig.0] For a 
periodic motion this curve is a closed loop and this is the 
case so for small d% = 0.02. However for large d% = 0.06 
the motion looks very complicated. The temporal evolu¬ 
tion of the spectrum is also not periodic for ds = 0.06, 
see panels (b) and (c). 

This behaviour can be explained by the fact that the 
pulsations of the synchrotron radiation causes pulsating 
recoil changing the intensity and velocity of the soliton. 
The variations of the soliton parameters in their turn af¬ 
fect the instantaneous frequencies of the radiation. For 
large d% the interplay of these effects leads to a complex 
dynamics of the solitons. Possibly this explains the com¬ 
plicated dynamics of cavity solitons discussed in [25| . 

Now we summarize the main results of the paper. It 
is shown that oscillating solitons can emit resonant radi¬ 
ation analogous to the synchrotron and cyclotron radia¬ 
tions of moving charges. The resonance condition for the 
radiation is derived and it is demonstrated that the pre¬ 
dicted resonances match the observed spectral lines very 
precisely. The structure of the emitted radiation and the 
effect of the recoil of the radiation on the solitons are 
discussed. The reported results are of general physical 
nature and thus can be observed not only in fiber cav¬ 
ities but in a wide class of physical systems supporting 
oscillating solitons. 
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